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Abstract 

The influence of the admixture of the state onto collective spin oscillations and neutrino 
emission processes in the triplet superfluid neutron liquid is studied in the BCS approximation. 



The eigen mode of spin oscillations with w ~ y^58/35A is predicted to exist in the triplet superfluid 
neutron condensate besides the already known mode uj ~ A/\/5- Excitation of the high-frequency 
spin oscillations in the condensate occurs through the tensor interactions between quasiparticles. 

Neutrino energy losses through neutral weak currents are found to consist of three separate con- 
tributions caused by a recombination of broken Cooper pairs and by weak decays of the collective 
modes of spin oscillations. Neutrino decays of the low-frequency spin waves can play an important 
role in the cooling scenario of neutron stars. Weak decays of the high-frequency oscillations that 
occur only if the tensor forces are taken into account in the pairing interactions does not mod- 
ify substantially the total energy losses. Simple expressions are suggested for the total neutrino 
emissivity. 



1 



I. INTRODUCTION 



It is considered well established that the inner core of neutron stars contains a condensate 
of superfiuid neutrons below the critical condensation temperature Tc. A superfiuid energy 
gap A arising in the quasiparticle spectrum suppresses most of the neutrino emission mecha- 
nisms in the volume of the star [I], especially when the temperature falls substantially below 
the critical value. In this case the number of broken Cooper pairs rapidly decreases and leads 
to a strong quenching of the neutrino emission caused by pair breaking and formation (PBF) 
processes, which is considered as the most efficient cooling mechanism of superfiuid neutron 



cores 



-5|. According to this scenario, at temperatures T ~ O.lTc, the neutron star enters 
the epoch of a surface cooling. 

The neutron superfluidity in the inner core of neutron stars is believed to arise owing to 
pairing of fermions into a triplet state. It is natural to expect the existence of low-frequency 
collective modes associated with spin fluctuations of such a condensate. Previously spin 
modes have been thoroughly studied in the p-wave superfiuid liquid ^He with a central 
interaction between quasiparticles 6|-l9| . These results cannot be applied without revision 
to the triplet superfiuid condensate of neutrons, where the pairing occurs mostly owing to 
a short-range negative spin-orbit force of the interaction in the channel of j = 2. 
In a series of papers 



10 



- |l2| we have investigated the collective spin oscillations in the 



superfiuid neutron liquid which can be formed because of the strong attractive spin- 
orbit interaction between neutrons at high densities. Spin waves with the excitation energy 
oj = A/\/5 were predicted to exist in such a superfiuid condensate and it has been shown that 
the spin-wave decay (SWD) through neutral weak currents leads to a substantial neutrino 
emission at the lowest temperatures T <^ Tc, when all other mechanisms of the neutrino 
energy losses are killed by the superfluidity. 

In this paper the problem is considered for the case of ^P2—^F2 pairing. We consider 
the spin eigenmodes of the superfluid condensate in the case of pairing owing to spin-orbit 
and tensor forces. The neutrino emission owing to PBF and SWD processes is calculated. 
The calculations are made within the BCS approximation by assuming a projection of the 
total angular momentum of the bound pairs rrij = as the preferable one at supernuclear 
densities. 

The paper is organized as follows. Section II contains some preliminary notes and outlines 



some of the important properties of the Green's functions and the one-loop integrals used 
below. We recollect the gap equations for the case of spin-orbit and tensor pairing forces. 
In Sec. Ill we discuss the renormalizations which transform the standard gap equations to 
a very simple form valid near the Fermi surface. In Sec. IV we derive, in the BCS approx- 
imation, the equations for anomalous three-point vertices responsible for the interaction of 
the neutron superfluid liquid with an external axial-vector field. In Sec. V we apply the 
angle average approximation to make the equations solvable analytically. In Sec. VI we 
analyze the poles of anomalous vertices to derive the dispersion of spin-density oscillations 
in the condensate. In Sec. VII we calculate the linear response of the superfluid neutron 
liquid onto an external axial-vector field. In Sec. Vlll we derive neutrino losses caused by 
the recombination of broken Cooper pairs and by the decay of spin waves. In Sec. IX, we 
evaluate neutrino energy losses in the ^P2— ^-^2 superfluid neutron liquid. Section X contains 
a short summary of our findings and the conclusion. 

Throughout this paper, we use the standard model of weak interactions, the system of 
units h — c—1 and the Boltzmann constant ks — 1- 

II. GENERAL APPROACH AND NOTATION 
A. Green functions and loop integrals 

The order parameter, D = Dap, arising due to triplet pairing of quasiparticlcs, represents 
a 2 X 2 symmetric matrix in spin space, {a, (3 ^t^J-)- The spin-orbit interaction among 
quasiparticlcs is known to dominate in the nucleon matter of a high density. Therefore it is 
conventional to represent the triplet order parameter of the system D — X^j/mj ^jimj^a^^^^ 
as a superposition of standard spin-angle functions of the total angular momentum (j, rrij), 



The angular dependence of the order parameter is represented by the unit vector n — p/p 
which defines the polar angles {9, cp) on the Fermi surface. 

For our calculations it is more convenient to use vector notation that involves a set of 
mutually orthogonal complex vectors hji„i. (n) defined as 




(1) 
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where <t = ((3"i, 0"2, (T3) are Pauh spin matrices, and g = i(72- The vectors hjim^ obey the 
normahzation condition 

Using the vector notation the order parameter is D (n) = Ah&g, where the vector b in 
spin space is defined by the relation 

A (p) b (n) = 5^ Ajim, (p) b,/^^. (n) . (4) 

jlmj 
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14j , b (n) is a real 



Because the ground state order parameter is to be a unitary triplet 
vector which we normalize by the condition 

^P(„) = l. (5) 

Making use of the adopted graphical notation for the ordinary and anomalous propa- 
gators, G = — »— , G~{p) = — • — , F^^^ = » « , and F^^^ = « > , it is convenient to 
employ the Matsubara calculation technique for the system in thermal equilibrium. Then 



the analytic form of the propagators is as follows 15|, llGj 



p) = G(p^,p)(5„^, G iPr„p) = G (p^,p)(5«/3, 
(p„ p)=F (p„ p) hag, F^^) (p^, p) = F (p„ p) gah, (6) 

where the scalar Green's functions are of the form G~ {prj, p) = G {—Pr], — p) and 

Here = ivr (2?7 + 1) T with = 0, ±1, ±2... is the Matsubara's fermion frequency, and the 
quasiparticle energy is given by 

El = el + A'b' (n) , (8) 

where Ep is the single-particle spectrum of the normal Fermi liquid. Near the Fermi surface 
one has 

The effective mass of a neutron quasiparticle is defined as M* = pp/vp, where <^ 1 is 
the Fermi velocity of the nonrelativistic neutrons. In the absence of external fields, the gap 
amplitude A (T) is real. 



The following notation is used below. We denote as Lx,x iy, q; p) the analytical contin- 
uation of the Matsubara sums: 

Lxx' (^w«,p+|;p-^) =T^X {p^ + u^,p+^X' (^p^,p-|) , (10) 

v 

where X, X' G G,F,G~, and = 2mKT with k = 0, ±1,±2.... We divide the integration 
over the momentum space into integration over the solid angle and over the energy according 
to 

[ d^p f dnl f°° 

J = ' J L"''-' 

where g = ppM* /n'^ is the density of states near the Fermi surface in the normal state, and 
operate with integrals over the quasiparticle energy: 

Xxx' (w, n, q; T) = ^ j dSpLxx' (w, p+|, P-^) • (12) 

These are functions of u, q and the direction of a quasiparticle momentum n. 

The loop integrals f[T2]) possess the following properties which can be verified by a straight- 
forward calculation: 

^g-g = ^gg- , ^gf = —^fg , '^g- f = —'^fg-, (13) 
Ig-f + Ifg = ^1ff, (14) 

^G-F — 1-FG = —^^FF- (15) 



For arbitrary u, q, T one can obtain also 

2A2 



Igg- + b-'lFF = A+ ^ ^ Iff, (16) 



where v is a vector with the magnitude of the Fermi velocity vp and the direction of n, and 

A (n) ^ [Igg- (n) + P (n) I^f (n)] . (17) 

B. Gap equation 

The block of the interaction diagrams irreducible in the channel of two quasiparticles, 
^ai3,'yS, is usually generated by the expansion over spin-angle functions ([1]). The spin-orbit 
interaction among quasiparticles is known to dominate at high densities. In this case the 



spin s and orbital momentum 1 of the pair cease to be conserved separately. Thus the 
complete list of channels participating in the triplet-spin P-wave pairing includes the pair 
states with j = 0, 1,2, and \mj\ < j. The pairing occurs in the state with j = 2 because 
the attractive interaction in this channel strongly dominates. The tensor components of 
the neutron-neutron interaction are known also to exert some influence on pair formation 
in dense neutron matter, favoring the condensation of pairs in the ^P2+^F2 state, but the 
contributions from — t-^Pq ot — t-'^Pi transitions are deemed to be unimportant. Hence 
we take the approximation to neglect the j = 0, 1 coupling throughout this paper. From 
now on we omit the suffix j everywhere by assuming that the pairing occurs into the state 
with j = 2. Then, in the vector notation, the pairing interaction is of the form 



^r,^,,, (p, pO = J2 (-1)"^ V,, {p,p') {himM)^9)^p (^^b*^^.(n')) (18) 

111.. • 



l-l' 

(-1] 

Vlmj 

where the pairing matrix elements Vw {p,p') with /, /' = j ± 1 = 1, 3 are the corresponding 
interaction amplitudes. 

The ground-state problem is normally treated in terms of the set of equations for the 



coupled partial- wave amplitudes A^^, 17H22|. Making use of the identity 



J^tanh% ^tV — (19) 
2E^ 2T Y^' + ^P 

one can obtain the standard set of equations for the triplet partial amplitudes Aimj in the 
form 



l'=1.3 ^ 



X A (pO ( h*,^^ (nOb(n')T ^ ) , (20) 

Here and in what follows we use the angle brackets to denote angle averages, 

(...) = — [ dn.... 

Notice that 

= G (p„ p) G- (p„ p) + ¥F ip„ p) F (p„ p) , (21) 

Pri + Fp 



and the gap equation (120!) can be identically written as 

;br,„^,(n')b(n') [Lgg- + b'L^p]^^,.) • (22) 



X 



III. RENORMALIZATIONS 



Both the gap equation fl22p and the vertex equation 0351) involve integrations over the 
regions far from the Fermi surface while we are interested in the processes occuring in a 
vicinity of the Fermi sphere. To get rid of the integration over the far regions we renormalize 



the interaction as suggested in Refs. [2J, |25|]: We define 

I" 

= Vw ip,p') ^V,. {p,p") (/) V,n' ip",p') , 

I" ^ ^ 

where the loop L^J^} {p") is evaluated in the normal (nonsuperfluid) state. Using the identity 



A (p) f dn 



bL.b = 1 



v.. {p,p') = Vw {p,p') ^V,, ip,p") 

1// ^ 



^Ira, {P) J 4:71 ^ 

one can recast the above as 

I'' 

n ' 

Then it can be shown [ll] that we may everywhere substitute V/^ for Vw provided 
that at the same time, we understand by the Lqq- element, the subtracted quantity 
Lgg-—L^q^_^ {p") [L^qq- is to be evaluated for a; = 0, q = in all cases]. The gap equation 
fl22l) becomes of the form 



Aim, = -5^^'-''V,,A(K,„^(n)b(n)A(n)) (23) 

I' 

which is valid in the narrow vicinity of the Fermi surface, where the smooth functions 
Aim, (p), {p,p'), and A{p) may be replaced with constants A (p) ~ A{pp) = A, etc. 
The function ffTTI) is now to be understood as 



A(n) 



IgG- — ^GG~ + ^ -^^P 
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(24) 

J a;=0,q=0 



It can be found explicitly by performing the Matsubara summation: 

.^^'^r,f 1 , + A262 1 e\ 

A(n) = - de \ , tanh tanh — . (25) 

^ ^ 2 + A252 2T e 2Tj ^ ' 

The renormalization of Eq. fl35|) also reduces to the replacements V^^' — V;;/. The function 
A (n) should be replaced by the expression (123]) . 



IV. VERTEX EQUATIONS 

We are interested in the linear medium response onto an external axial-vector field. The 
field interaction with a superfluid should be described with the aid of two ordinary and 
two anomalous three-point effective vertices. In the BCS approximation, the ordinary axial- 
vector vertices of a nonrelativistic particle and a hole are to be taken as a and <T"^, re- 
spectively. The anomalous effective vertices, T^^-* (n;a;, q) and T'-^^ (n;a;, q) are given by the 
infinite sums of the diagrams taking account of the pairing interaction in the ladder approx- 
imation [23]. These 2x2 vector matrices are to satisfy the Dyson's equations symbolically 
depicted by graphs in Fig. [H 



y^^>. y^^y 



FIG. 1. Dyson's equations for the anomalous vertices. The ordinary vertices are shown by small 
filled circles. The shaded rectangle represents the pairing interaction. 



The analytic form of the diagrams in Fig. 1 is derived in Ref . [10| . We are interested in the 
neutrino energy losses through neutral weak currents. In the case of a nonrelativistic medium 
the relevant input for this calculation are the effective weak vertices at zero momentum 
transfer. This substantially simplifies the problem. After some algebraic manipulations the 
BCS equations for anomalous vertices at q = can be found in the following form (we omit 
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for brevity the dependence of functions on u): 



(n) = J2^^irr.,H9j2^u'^ (^GG-Tr [g (^b*„J t^^ 

Inij I' 

-Xp^Tr (^b) ^t(2) (^b)] - ^X^^2z (b^^^ xb 



(26) 



t(2) (n) = 5^^<TbL^,(n) 5^ V.,- (Xc-cTr [(^b,„^,) 



.(2) 



-X^i^Tr [i&hi,^^] (o-b) t(i)^ (a-b)J - ■^Xi.^2z (b^^^^.xb) 



(27) 



Inspection of the equations reveals that the anomalous axial- vector vertices can be found in 
the following form 

t«(n,a;) = 5^B« (a;)(^b,^J^, (28) 



Irrii 



Irrii 



(29) 



Insertion of these expressions into Eqs. fl26|) and f l27|) makes it possible to obtain the following 
equations 



T3 + 



l"m.'. 



-2 J2 { (bbL^.) (bb,,„, ) Iff) B+„, + | (x^^^ (b;,„^ xb 



(30) 



B 



Im-i 



J2^u'J2[{i^oG--b'iFF) (br,„,b 



l"m'. 



l"mr 



+2 bb,„ J bb,w X^^ 



B 



l"m'. ' 



where the new unknown vector functions are defined as 



B 



(±) 



Im 



J 2 



B 



(2) 



(31) 



(32) 



The uniform Eqs. ( 1311) have nontrivial solutions if the determinant of the system equals 
zero. This condition could be considered as the dispersion equation for the eigen modes of 
oscillations. However, the solutions obtained in this way would be spurious, because the 
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physical solution must develop a pole in the vertex function when the frequency approaches 
the eigenvalue. In the case of uniform equations the solution remains finite at the resonant 
frequency. Therefore only trivial solutions B^^, = are physically meaningful. We then 
obtain 



Irrii Irui 



mj -0(2) 



(33) 



From now on we omit "plus" in the notation by assuming B^^j = B^,. 

By making use of Eq. ( TT6l) and denoting Xiri? (n, q = 0) = Xq (n,a;; T), where 



Xo (n,w) = / — ^ 



AEl - (w + iO) 



2 2T' 



(34) 



one can obtain Eq. fl30|) in the form 



B 



2A2 



2 bb,„ J bb 



+z-(Xo(n) (K,„^.xb 



(35) 



This equation is to be solved together with the gap equation 



V. ANGLE AVERAGE APPROXIMATION 

The angle dependence of the functions A{ri) and Xpp{n) arises owing to anisotropy 
of the square of the energy gap (p') 6^ (n) entering the energy of a quasiparticle. The 
consideration can be substantially simplified using the angle average approximation, that 
is, replacing the anisotropic energy gap with its angle average, A^6^ — )■ (|A^6^) = A^. 



In Refs. 



26 



28|, it has been shown that the angle average approximation is an excellent 
approximation to the true solution, as long as one is only interested in the average value of 
the gap near the Fermi surface, and not the angular dependence of the gap functions. After 
this replacement the angle integration becomes trivial. Making use of the orthogonality 
relations ([3]) after the renormalizations from Eq. (|35|) we get a set of linear equations for 
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each value of rrij 



- 2 ^i"m', ((bbL,) (bb,,^,))x3v +z|Xav (br,„,xb)} , (36) 



where 



and 



1 p /I E 1 F . 

^av = - / I ;^ tanh ^ - - tanh — ) , (37) 



Xav = / ^ ^ tanh — (38) 







with E = Ve^ + A2. 

The gap equation (l23l) becomes of the form 

A,„^, = - ^ ^'-''Vi^ A^^^A^v (39) 

I' 

It is convenient to write the Eqs. (136|) and ( 139|) as two matrix equations. The correspond- 
ing vertex equation is 

^av ~l~ 2A^ -^av^ Bx^ij 

2Xav E.™; ( (bbLj (bb,„, ) ) B,„, ^ ^ . Xav (bU^. xb/ I ^^^^ 



/ Bim^ 1 










2Xav E.n; ( (bbL^ ) (bb,^ ) ) B,„, y ^ \ Xav (b^„^ X b 

and the gap equation becomes of the form 





(41) 



In obtaining the equations the fact is used that the interaction matrix is symmetric on the 
Fermi surface, V31 = V13. 

The interaction matrix which enters Eqs. fHOl) and (HTl) can be diagonalized by unitary 
transformations V' = UVlf^ with U being the unitary matrix 



U={U-^y = - ^1 ^'^ 1, (42) 




where = ^{V^^ - Vn)' + 4V?3 ± (V33 - Vn). 
One has UVW = diag W+) with 



= ^ ( V33 + Vn ± \/(V33 - Vn)' + 4Vf3 ) . (43) 



Applying the unitary transformation U to the gap equations f l4T|) yields two coupled equa- 
tions: 

^/V^Al^^ + ^/KA3rn, = -W^ (v^Ai„^ + ^/vlAsm,) (44) 

yV^Ai^^ - ^/V^Asrn, = -W+ (^^/KAi^^ - y^Asrn,) ^av (45) 

In Eq. fj40|) . the interaction matrix can be also diagonalized by the unitary transformation 
Further simplification is possible owing to the fact that by virtue of Eqs. (jH]) and 
the coupling constants W± can be removed out of the equations. This results in the 

set of equations 

- E (v^((bbL,) (bb,„;)) + V^((bbL^.) (bb.„;))) B,„, 

= {^'-^ X b) + (b;^^, X b)) , (46) 



J2 (v^((bbL^) (bb,„,)) - v^((bb;„^) (bb,„,))) B,^, 



= ~Ya' (^^-^ X b) - (bL, X b)) . (47) 

Equations ( 146|) . and (147|) represent a closed set of equations for determining the vertex 
amplitudes B^^ (i^)- 



VI. EIGENMODES OF SPIN OSCILLATIONS 

For further progress we have to define the ground state of the condensate which is specified 
by the vector b(n). We focus on the condensation with vrij = which is conventionally 
considered as the preferable one in the bulk matter of neutron stars 

Q, a, a. 

In this 

case only the partial gap amplitudes with / = 1,3 and rrij = contribute. To simplify the 
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notation we denote them as Aio = Ai and A30 = A3. Taking into account spin-orbit and 
tensor interactions the ground state of such a triplet condensate is given by the vector b (n) 
of the form 

bfn) = ^bio(n) + ^b3o(n), (48) 



where A^ = A^ + A3. For this particular form of the ground state the following relations 
can be verified by a straightforward calculation 



(bb;„^ ) ( bb^,. 



5m,m'. (bb^^J (bbp^.) . 



We denote 



and 



bb* (bb,.^J 



2A 



(49) 



(50) 



From Eqs. ( I46ll and ( l47ll one can obtain five sets of linear equations corresponding to 
different values of nij = 0, ±1, ±2. 

For rrij = we obtain a set of two equations: 



+ 



Bio 
B30 = 0, 



(51) 



l^(fi2-/3o''') + v^A 



Bio 

B30 = 0, 



The sets of equations for m,- = ±2 are of the form 



(52) 



+ 



ir(fi2_^3,3) _y^^l,3 



Bl,±2 
B3,±2 = 0, 



(53) 



B 



3, ±2 



0. 



(54) 
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For each of the values of rrij = ±1 we get a set of two equations: 



1,3 

±1 



Bi,±i 



±1 



V- (b* ±1 X b 



(55) 



3,1 
±1 



Bi 



±1 



B 



3,±1 



-Qi 



^-(bUiXb>-v^(b*,±iXb> . 



(56) 



By the same reason as in the case of Eqs. (|3T|) the uniform equations for m,- = and for 



±2 have only trivial physical solutions Bio = B30 = and Bi^-1-2 = B3^±2 = 0. The 



solutions to Eqs. f l55|) and fl56|) are 

-in 



Bi,±i 

B3,±l 



X±i (fi) 

-in 
x±i m 



where 



[{n' - {hl^, X b> + (b^i X b>] 
[{Q' - {hl^, X b> + (b^i X b>] 
1 /3 Ai 



(bUiXb> = -^v^^ (^,±1,0) 



(57) 
(58) 

(59) 
(60) 



and 



(61) 



Notice that the interaction parameters V± drop out of the final result in Eqs. fl57|l and fl58l) . 

Poles of the effective vertices at x±i (^) = signal the existence of collective modes. 

Frequencies of the eigenoscillations are independent of the sign of rrij. For each of the values 

of rrij = ±1 we obtain two eigenmodes: 

1 
2 

(62) 



Straightforward calculations give 



I A'W 2 ^As 58A2 

--^'^-r^Al + W 



20 A2 



(63) 
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It is interesting to compare the resonant frequencies fl62|) with the spin oscillation fre- 
quency ojs = A/-\/5 obtained in Refs. 10|, luj in a simple model restricted to excitations of 
the condensate with / = 1. By taking A3 = and Ai = A in the matrix elements we find 
4/3^f ^ — By neglecting the small term 4/3^f under the root in Eq. 



we obtain 



ui (A3 = 0) = 2A/3]ll (A3 = 0) = ^A ~ 0.45A, (66) 

V5 



and 



0J2 (A3 = 0) = 2A/3if (A3 = 0) = y^A ~ 1. 29 A. (67) 

We see that the extending of the decomposition scheme of the excited states with the total 
angular momentum j = 2 up to / = 1, 3 leads to a very small frequency shift of the known 
mode, cOs = A/v^, but opens the new additional mode of collective spin oscillations. Inclu- 
sion of the tensor interaction implies A3 7^ and A^ = A^ + A|. In this case from Eq. (162|) 
we obtain two twofold {rrij = ±1) frequency: 



2 ,2 A3 1 ^A3 23 A2 
^-^nii-3^A^ + 21Af 



,15 25 /-A3 485Ai 370 /-Ai 55A^, 
28 ~ W^A; + 147 Af " 441^^ Af + 63 Af ' ' ^^^^ 



^Vl4 3'^Vi + 21A? 



/15 25 /-A3 485 AI 370 r-A^ 55 A^ 
28 ~ 49^ A^ 147 Af " 441^ Af + 63 Af 



(69) 



According to ca' 
(see, e.g., Ref. 
with a; = cji ~ 



culations of different authors, at the Fermi surface one has A3 ^ 0.17Ai 
isl). In this case our theoretical analysis predicts two degenerate modes 
0.42 A and two degenerate modes with u = U2 = 1. 19 A. 
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VII. ANOMALOUS VERTICES AND POLARIZATION FUNCTIONS 



Making use of Eqs. (12^]). and (15^ we find 



T^^) (n, cj) = ^ [Bim^. (a-him^) g + 63^^. (a-hsm,) g] 



(70) 



(71) 



where the functions Bi^ (w) and Bsm (w) are given in Eqs. fl57|) and f l58|l . 



The general expression of the axial polarization tensor in 
been already discussed before. It can be obtained in the form 



;he BCS approximation has 



lomJl 



gVLl^y ( Tr 



Tr 



0", 



(^b) gT^ 



(2) 



(72) 



where i,j = 1, 2, 3 and the anomalous axial- vector vertices T*^^'^^ are given by Eqs. ( !70|) and 
( I7T|) . We omit for brevity the dependence on n and u. Calculation of the traces results in 
the expression 



Irrij X 



(73) 



Inserting the vectors B^ ,„ , given by Eqs. ( 157|) . and ( 158|) . we write the function ( 16T|) . as 



(74) 



where rrij = ±1, and 



4A2 ' 



^2 

4A2' 



(75) 
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In this way we obtain 



Q2 



X 



E 



+ p'J^{h,m,y<hy (bL, xb 

+ <b3^.xby(bL, xb 

+/5^;<b3n.,xby(bL^. xSyj. (76) 

The poles location on the complex plane of fl is chosen so that to obtain the retarded 
polarization function. 

Summation over nij — ±1 can be done using the fact that the parameters entering 
this equation are independent of the sign of mj, and = P^l- Then a simple calculation 
gives 

J2 (bzm, xb>' (b;,^. X b)' = {Sij - SisSja) 



mi=±l 



(77) 



with 



and 



3A2 lAi IA1A3 

All - , A33 - , Ai3 - A31 - 



1 A^ 

We finally obtain the expression 



12 A^ 



6A^ 



(78) 



(79) 



X 



All 



u 

4A2 



+A 



a;2 - ^2 + ^0) (a;2 - a;| + iO) 
2 



33 



1 + 2Ai3^| 



4A2 



1,3 



(80) 



Below we use the retarded polarization tensor for a calculation of the neutrino emissivity 
of a nonrelativistic superfluid matter. In this calculation one can neglect the temporal and 
mixed components of the polarization tensor occurring as small relativistic corrections. 
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VIII. NEUTRINO ENERGY LOSSES 



We examine the neutrino energy losses in the standard model of weak interactions. Then 
after integration over the phase volume of freely escaping neutrinos and antineutrinos the 
total energy which is emitted per unit volume and time can be obtained in the form (see 



details, e.g., in Ref. |29| ) 

e = -^f^ f dco I rf^^^lly^ Imn- (c, q) {k,K - A: V) , (§1) 

where Gp is the Fermi coupling constant, Ca = 1-26 is the axial-vector weak coupling 
constant of neutrons, AC = 3 is the number of neutrino flavors, O (x) is the Heaviside step 
function, and = {u, q) is the total energy and momentum of the freely escaping neutrino 
pair (yU, u = 0,1, 2, 3). 

In Eq. ( ISTl) . we have neglected the neutrino emission in the vector channel, which is 
strongly suppressed owing to conservation of the vector current js^, [l^. Therefore the 
energy losses are connected to the imaginary part of the retarded polarization tensor in the 
axial channel, Imn^*^ ~ 6^'^6'^^ ImlV^. The latter is caused by the PBF processes and by 
the SWDs. These processes operate in different kinematical domains, so that the imaginary 



part of the polarization tensor consists of two clearly distinguishable contributions, Im 11^ 

•■PBF ^ -^^^^-^-^SWD; 



Imllppip + Imllowp,, which we now consider. 



A. PBF channel 

The imaginary part of Xgy, which arises from the poles of the integrand in Eq. ( 138|) at 
|a;| = 2E is given by 

Iml.. > 0) = - t^h (-) . (82) 

With the aid of this expression we find 

A20 (tj2 _ 4A2) / a; \ 
Imnr(.)^-2., ^X-_,^.^ anh(-) 

X <^ Sij - {bibj) - {6ij - 6i363j) -- 



2^ 



X 



4 A2 (a;2 _ ^2) (^2 _ ^ 

uP_ _ 3,3\ 4 A| _ , A _ 4 A3 

4A2 ^1 J 9AI V4A2 ^1 J 3A/1 



(83a) 
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Inserting the imaginary part of the polarization tensor into Eq. (IHTl) . we calculate the 
contraction of Imllpgp with the symmetric tensor — k'^g^v This gives 

1 1 



epBF 



1 /"°° 



967r6 



X 



, , , tanh ( — z- ) 

2A _ 4^2 exp (^) - 1 \4T ) 

18 A2 . , 



q<ijJ 

3 
' 4 



6 A2 


rA^12A| 
A 7 A? 


A? 


4A2u;2 



7 A2 



A2 (a;2 _ ^2) (^2 _ ^2) 



X 



UJ 



Pi 



1,1 



4 As .1,3' 
3A;^^ 



(84) 



4A2 y ' 9A2 V4A2 

Integration over d^q can be done in cylindrical frame, where qi = q± cos $, 52 = (1± sin 
and qs = qz- This results in the neutrino energy losses of the form 



epBF 



4 3 A2 



+ expz) 



X 



5 5A2(x2+?/2(l_fi2))(3;2+^2(l_^2)) 

,2 



I 2/1 o3,3\ I 4A3 . 2 , 2/1 pl,l\\ 4 A3 2ol,3 

+ 1/ ) + g^(a; +y [l - /3{ )) - -—y P{ 



(85) 



where z = ^/x^~+~y^, y = A (T) /T, and ^1^2 are defined in Eq. fl75|) . In obtaining Eq. fl85|) 
the change is used w = 2T ^ + jT^. 

For a practical usage from Eq. f l85|) . we find 



epBF = 5. 85 X 10 
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erg 
crrr^s 



(86) 



where M is the bare nucleon mass; Tg = T/10 K, and 

v4 



-^PBF (y) = Z/M dx- "2 

Jo (1 + exp z) 

.2 



X M-3 



A? 



a;2 + ?/2 



A2 (s2 + ?/2 (1 _ f]2)) ^^2 J^yi{\- ^2)) 



(87) 



In the limit A3 = 0, the neutrino energy losses, as given by Eq. fISS]) reproduce the result 



obtained in Ref. 



12| for the one-component phase vnj = 0. It is necessary to notice that Eq. 



E7} obtained in the angle-average approximation is much simpler for numerical evaluations 



than the "exact" expression which contains additionally the angle integration 
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ll|. To avoid 



possible misunderstanding we stress that the gap amplitude A (T) in Eq. (157]) is \/2 times 



larger than the gap amplitude ^ykl used in Ref. 
An = A6 (n) is written in the form An = Ay^^Vl + 3 cos^ 6 
words, (A^) = A^ 



31| , where the same anisotropic gap 
Ayi^-LV^ 6 (n). In other 



'^'-^YKL- 



B. SWD channel 

In the frequency domain < w < 2A, the imaginary part of the weak polarization tensor 
( IHOj) arises from the poles of the denominator aX u = uji and uj = UJ2 and consists of two 
terms 



Imn^WD > 0) 



A2 [nl-nl] 



4A^ 



4 A 



3 Ai 



(88) 



According to Eqs. (ETD at A3 = one has 0.2 = /^If and ImH^w^ {u U2) = 0. In other 
words, the high-frequency spin oscillations can not be excited if the tensor interactions 
between the pairing particles are not taken into account. 

Inserting Eq. f l88|) into Eq. ( 18T|) and performing trivial calculations, we find two con- 
tributions to the neutrino energy losses. The first contribution is caused by the decay into 
neutrino pairs of the lowest mode of spin oscillations at co = ui. 



"SWD 



3207r5 



3,3 



A2 Qj- 
4A2 



3a/^ +9 a? ^ 



X 



col 



exp if) - 1 



de 



A' 



, E 
„ , - tanh — . 
E E^-ujI/4: 2T 



(89) 



According to this equation, in the case of A3 — )• 0, Ai — )• A, and ooi — t- A/ 75, the energy 
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losses are twice less than that found in Ref, 



e« D (A3 = 0) = ^^GlClKppM 



ool de A" ^ E 



X / ^tanh— . (90) 

exp(f ) -1 io E E^- col/A 2T ^ ' 



We use the opportunity to point out the error in Eqs. (79) and (81) of Ref. ([12]), where 
the factor of 1/2 is lost. 

The second contribution originates from weak decays of the second (higher) mode at 
UJ = U2: 



cul de A2 , , E 



tanh— . (91) 



exp(f) - Uo E E'^-ujI/A 2T' 

Because excitation of the high-frequency spin oscillations occurs through the tensor in- 
teractions, the contribution of the second mode vanishes if the tensor forces are switched off 
(i.e. when A3 = 0) . 

The expressions (l89|) and (lOTl) can be written in the traditional form 



e^') -1 

^SWD ~ ^ 



-J'swDi^uy) (92) 



g(2) 



Ql — Ql exp (2yi7i) — 1 ' cm^s' 

'^nFn! exp(2yn,)-i -^^'"°'"^-^' 

where /Sj'^ , given by Eqs. (1631) -(l65l). are functions of the gap components Ai and A3; 
y = A{T) /T, and 



00 



du 1 y 



^swD {n, y) = y' -tTtt 2^1 02 ^^^^ f V^^- (94) 
Jo Vu^ + 1 u-^ + 1 - iV 2 
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26 



rrij =0 




iog,.T [K] 



FIG. 2. Temperature dependence of the neutrino emissivity owing to recombination of Cooper 
pairs (PBF) and owing to decay of spin waves (SWDi, the emissivity of the lower mode. SWD2, 
the emissivity of the upper mode) at A3/A1 = 0.17 and pp = 2.1 fm~^. 

IX. EFFICIENCY OF THE NEUTRINO EMISSION 

In general, the temperature dependence of the gap amplitudes Ai and A3 is to be found 
with the aid of the gap equations. For simple estimates we take the approximation that the 
ratio A3/A1 remains constant when the temperature varies, and the temperature dependence 
of the gap is given by the function y = A (T) /T. This function is well investigated for a 



3 



P2 pairing. Since the tensor contribution can be considered as a perturbation 2l|, in a 



zero approximation, we can use, for example, the simple fit to Aykl (T) /T = vb (t), as 
suggested in Ref. |3l|, where r = T/T^. Taking into account that, in Ref. 3l|, the gap 
amplitude Aykl (T) is defined by the relation A^ = Ayj^j^ (1 + 3 cos^ 9), while our definition 
is Al (A3 = 0) = |A2 (1 + 3 cos2 9), we obtain y (r) = v^v^ (r). 

In Fig. |2]we show the neutrino emissivity e caused by the PBF processes and by the decay 
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of the lowest mode (SWDi) and the higher mode (SWD2) of spin oscillations. The temper- 
ature dependence of the emissivity is evaluated at pp = 2.1 fm~^ assuming A3/A1 = 0.17. 
We set the effective nucleon masses M* = 0.7 M; the critical temperature for neutron pairing 
is chosen to be = 3 x lO^K and = 3 x lO^K. 

One can see that the decay of the low-frequency spin waves into neutrino pairs (SWDi) is 
very effective at low temperatures, when other known mechanisms of neutrino energy losses 
in the bulk neutron matter are strongly suppressed by the superfluidity. As discussed in Ref. 
[12 1 the neutrino emission caused by the decay of the low-frequency spin- waves can dominate 
the 7 radiation within a wide range of low temperatures, which was considered before as the 
photon-cooling era. A simple estimate has shown that the decays of spin waves can modify 



the cooling trajectory of neutron stars (see Fig. 5 in Ref. 12|). 

Weak decays of the high-frequency mode of spin oscillations occurs only if the tensor 
forces are taken into account in the pairing interaction, that is, if A3 7^ 0. Although the 
maximal neutrino emission caused by the SWD2 processes is as large as in the SWDi the 
neutrino energy losses from the decay of the upper mode decrease more rapidly along with 
lowering of the temperature. As a result the SWD2 contribution into the total energy losses 
is negligible in comparison with the sum of the PBF and SWDi contributions. We found 
that the latter can be excellently described by the expressions obtained in Ref. for the 
case of pairing with rrij = [However see note after Eq. fl90|) ]. We quote these simple 
expressions for references: 

.PBP = 5. 85 X 10- (^) (£) T^KClF,,, (y) (95) 



with 



and 



^4 



i^PBF (y) = / dx- — -2 , (96) 

Jo [l + expz) 

.SWD = 1. 37 X 10- ( !^) (If) TXCl-^^S!^ ^. (97) 

\MJ\McJ " *exp(!//v'5)-l cm3s ^ ' 



erg 
cm^s 



du y 



^° = / 7^-^ ^^^^ ^vvn. (98) 

Jo [u^ + 1) ^ 



where 



The accuracy of Eq. (|97|1 substantially increases if one takes into account that the spin 
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FIG. 3. Temperature dependence of the total bulk neutrino luminosity from a homogeneous super- 
fluid core owing to recombination of Cooper pairs (PBF) and owing to decay of spin waves (SWD) 
for pp = 2.1 fni^^. The effective mass is taken to be M* = 0.7M. Solid lines are calculated 
according to the exact expressions ([55|) . (^^. Dash lines are calculated by simplified formulas 
(1951) and (1971) . Volume of the triplet condensate is estimated as 7 x 10^*^ cm'^. The short-dashed 

n 

line is the energy losses per unit of time owing to the surface 7-radiation, as calculated in Ref. [4]. 



wave energy ui depends on the temperature. As found in Ref. [11] this dependence can be 
evaluated making use of the analytic fit which relates Qi to y at any y > 0: 



The maximum fit error is about 0.1%. 

To get an idea of the accuracy of the simplified expressions, in Fig. [3] we demonstrate the 
total neutrino energy losses caused by PBF and SWD neutrino emission from the superfiuid 
core of the volume 7 x 10^^ cm^. The total neutrino luminosity, as calculated by the exact 
Eqs. ( 186|) . ( 192|) . and (|93|) . is shown in comparison with the sum of the PBF and SWD 




wi _ 0.2172 - 0.005% + O.OllV + 0.0026?/=^ 
2A ~ 1 + 0.053% + 0.0710?/2 + 0.0175?/3 



24 



neutrino losses calculated with the aid of simple expressions given by Eqs. ( !95|) and (1971) . 
We display also the the luminosity of the surface photon radiation. The latter is taken as 
in Fig. 20 of Ref. fl- 



X. SUMMARY AND CONCLUSION 

According to modern theories the triplet-spin superfluid condensate in the inner core of 
neutron stars arises owing to pairing of neutrons caused by attractive spin-orbit and tensor 
forces and consists of the orbital contributions corresponding to / = 1,3. Assuming the 
projection of the total angular momentum rrij = , the superfluid energy gap A arising in 
the ^P2—^F2 state is = + A3, where the contribution A3 7^ is caused by the tensor 
interactions. We have studied the influence of the admixture of the state onto the 
collective spin oscillations and neutrino emission processes in the triplet superfluid neutron 
liquid. 

To evaluate the rate of neutrino energy losses out of the ^P2—^F2 superfluid neutron 
liquid we have calculated the anomalous three-point vertices responsible for the interaction 
of the superfluid liquid with an external axial-vector field. The BCS-like calculation has 
done in the angle average approximation. The polelike behavior of the vertices points out 
the existence of two twofold eigen modes of oscillations. The oscillation frequencies in 
terms of the gap components are given by Eqs. fl68|) and fl69l) . According to the obtained 
expressions the known low- frequency mode coi ^ A/-\/5 undergoes only a small frequency 
shift owing to the tensor interactions. The frequency of the new, upper mode U2 (A3 = 0) ~ 
a/58/35A decreases slightly along with increasing of the tensor contribution into the energy 
gap. We found that the excitation of the high-frequency spin oscillations is strongly quenched 
if the tensor interactions between the pairing particles are not taken into account, that is, 
if A3 = 0. According to calculations of different authors, at the Fermi surface one has 



A3 ~ O.I7A1 (see, e.g., Ref. [l8|). In this case our theoretical analysis predicts two twofold 
modes a; = wi ~ 0.42A and u = U2 = 1- 19A. 

We have derived the linear response of the superfluid liquid onto an external axial- vector 
field. At the time-like momentum transfer the imaginary part of this polarization function 
consists of three contributions originating from a recombination of broken Cooper pairs and 
from weak decays of the collective modes of spin oscillations. Accordingly, the neutrino en- 
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ergy losses through neutral weak currents consist of three contributions caused by the above 
decay processes. The neutrino energy losses owing to PBF, SWDi, and SWD2 processes are 
presented analytically by Eqs. ( M2\i . and (jHSD- 

Neutrino decays of the low-energy spin waves (SWDi) can play an important role in the 
cooling scenario of neutron stars. Previously we have demonstrated (see Fig. 5 in Ref. 12|) 
that the decays of spin waves with u = A/y/E can become the dominant cooling mechanism 
in a wide range of low temperatures and modify the cooling trajectory of neutron stars. 

Weak decays of the high-frequency mode (SWD2) occur only if the tensor forces are taken 
into account in the pairing interaction, that is, if A3 7^ 0. The maximal neutrino emission 
caused by the SWD2 processes is of the same order as in the SWDi, however the neutrino 
energy losses from the decay of the upper mode decrease more rapidly along with lowering of 
the temperature. As a result the SWD2 contribution into the total energy losses is negligible 
in comparison with the sum of the PBF and SWDi contributions. This fact makes it possible 
to neglect the SWD2 contribution and describe the neutrino energy losses from the ^P2—^F2 
superfiuid liquid by simple expressions given by Eqs. ( 195|) and ( 197|) . 
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